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The frequencies of the highest-frequency kilohertz QPOs recently discovered in some sixteen neutron stars in 
low-mass X-ray binary systems are most likely the orbital frequencies of gas in Keplerian orbit around these 
neutron stars. If so, these QPOs provide tight upper bounds on the masses and radii of these neutron stars and 
important new constraints on the equation of state of neutron star matter. If the frequency of a kilohertz QPO 
can be established as the orbital frequency of gas at the innermost stable circular orbit, this would confirm one of 
the key predictions of general relativity in the strong-gravity regime. If the spin frequency of the neutron star can 
also be determined, the frequency of the QPO would fix the mass of the neutron star for each assumed equation of 
state. Here we show how to derive mass and radius bounds, using the kilohertz QPOs, for nonrotating and slowly 
rotating stars, and discuss how these bounds are affected by rapid stellar rotation and radial radiation forces. We 
also describe observational results that would be strong evidence for the presence of an innermost stable circular 
orbit. No such strong evidence is present in current data, but future prospects are excellent. 



1. INTRODUCTION 

Determination of the equation of state of neu- 
tron stars has been an important goal of nuclear 
physics for more than two decades. Progress to- 
ward this goal can be made by establishing astro- 
physical constraints as well as by improving our 
understanding of nuclear forces. 

Many ways of deriving constraints from astro- 
physics have been explored. One of the best 
known is pulse timing of pulsars in binary sys- 
tems. Although binary pulsar timing has made 
possible stringent tests of general relativity (see, 
e.g., [16]), the w 1.4 masses derived from tim- 
ing (see [17]) are allowed by all equations of state 
based on realistic nuclear physics, and hence these 
observations have not eliminated any of the equa- 
tions of state currently being considered. The 
highest known neutron star spin frequency, the 
643 Hz frequency of PSR 1937-^21 [1], is also al- 
lowed by all equations of state currently under 
consideration. Radius estimates based on the en- 



ergy spectra of type I X-ray bursts and on obser- 
vations of thermal emission from the surfaces of 
neutron stars are more restrictive in principle but 
currently have large systematic uncertainties (see 
[6,11]). ^ 

The discovery of high-frequency brightness os- 
cillations from some sixteen neutron stars in low- 
mass X-ray binaries with the Rossi X-ray Timing 
Explorer (RXTE) holds great promise for pro- 
viding important new constraints. Oscillations 
are observed both in the persistent X-ray emis- 
sion and during type I X-ray bursts. The kilo- 
hertz quasi-periodic oscillations (QPOs) observed 
in the persistent emission have high amplitudes 
and relatively high coherences (see, e.g., [18]). 
A pair of kilohertz QPOs is commonly observed 
in a given source. Although the frequencies of 
these QPOs vary by up to a factor of ^2, the 
frequency separation Av between a pair of kilo- 
hertz QPOs appears to be constant in almost all 
cases. In both the sonic-point [10] and magne- 
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tospheric [15] beat-frequency interpretations, the 
higher frequency in a pair is the orbital frequency 
at the inner edge of the Kcplcrian flow, whereas 
the lower frequency is the beat of the stellar spin 
frequency with this frequency. Such high or- 
bital frequencies yield interesting bounds on the 
masses and radii of these neutron stars and inter- 
esting constraints on the equation of state of neu- 
tron star matter. In § 2 we describe how mass and 
radius bounds can be derived from the proper- 
ties of the kilohertz QPOs, and discuss how these 
bounds arc affected by rotation and by radial ra- 
diation forces. 

The bounds on the mass and radius and on the 
equation of state would become particularly re- 
strictive if the frequency of a kilohertz QPO can 
be securely established as the orbital frequency 
of gas at the innermost stable circular orbit. In- 
deed, if the QPO frequencies currently observed 
in some sources are the orbital frequency of gas 
at the innermost stable orbit, several currently vi- 
able equations of state are ruled out. Moreover, 
detection of the effects of the innermost stable 
orbit would, by itself, be a conflrmation of one 
of the key predictions of general relativity in the 
strong-gravity regime. Several authors have re- 
cently suggested that these effects have already 
been observed and have argued that the neutron 
stars in some kilohertz QPO sources therefore 
have masses close to 2.0 M©. In § 3 we discuss 
these suggestions, and describe observational re- 
sults that would be strong evidence for the pres- 
ence of an innermost stable orbit. We argue that 
the evidence cited to support detection of the in- 
nermost stable circular orbit was not compelling 
and that more recent observational results do not 
support these claims. 



2. CALCULATIONS 

Suppose that, as in the sonic-point model, the 
frequency fQpo2 of the higher-frequency QPO in 
a kilohertz pair is the orbital frequency of gas 
in a nearly circular Keplerian orbit around the 
neutron star and that the highest observed value 
of t'QP02 from a given star is i'qpo2- 



2.1. Nonrotating stars 

Assume first that the star is not rotating and is 

spherically symmetric. Then the exterior space- 
time is the Schwarzschild spacetime and the or- 
bital frequency (measured at infinity) of gas in a 
circular orbit at Boyer-Lindquist radius r around 

a star of mass M is 

i,0(M,r) = (l/27r)(GM/r3)i/2 . (i) 

Here and below the superscript indicates that 
the relation is that for a nonrotating star. Equa- 
tion (1) may be solved for the mass of the star as a 
function of the radius of the orbit with frequency 
'^QP02' with the result 

MO(iiorb, ^^QP02) = (47rVG)i?^,b('^QP02)' • (2) 

The mass of the star is related to the radius of 
the innermost stable orbit by the expression 

M°{R^,) = {c'/QG)R^, . (3) 

The radius of the star must be smaller than the 

radius i?orb of the gas with orbital frequency 
'^QP02) so the representative point of the star in 
the i?,M plane must lie to the left of the curve 
M'^(i?orb, f^QP02)- addition, in order to pro- 
duce a wave train with tens of oscillations, the 
gas producing the QPO must be outside the ra- 
dius R-ms of the innermost stable circular orbit, so 
the representative point must also lie below the 
intersection of M°(i?orb, '^qpo2) with the curve 
Af°(i?,„s). Figure la shows the allowed region of 
the R,M plane for J^Qpo2 = 1220 Hz . The maxi- 
mum allowed mass and radius are [10] 

il/C = 2.2 (1000 Hz /i.^po2) Mq (4) 
<ax = 19.5 (1000 Hz /z.^po2) km. (5) 

For example, the 1220 Hz QPO observed in the 

atoll source 4U 1636—536 would constrain the 
mass of this neutron star to be less than 1.8 
and the radius to be less than 16.0 km, if it were 
not rotating. Figure lb compares the mass-radius 
relations for nonrotating stars given by five equa- 
tions of state with the regions of the radius-mass 
plane allowed for three values of fc'Qpo2- 
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Figure 1. (a) Radius-mass plane, showing the region allowed for a nonrotating neutron star with 
^QP02 ~ 1220 Hz. (b) Comparison of the mass-radius relations for nonrotating neutron stars given 
by five representative equations of state with the regions of the mass-radius plane allowed for nonrotat- 
ing stars and three different Keplerian orbital frequencies. The light solid curves show the mass-radius 
relations given by equations of state A [12], FPS [7], UU [20], L [14], and M [13]. (c) Regions allowed 
for rotating neutron stars with various values of j and i^qpo2 = 1220 Hz, when first-order effects of the 
stellar spin are included, (d) Illustrative Keplerian QPO frequency given by fully general relativistic 
calculations of the gas dynamics and radiation transport in the sonic-point model [10]. 
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2.2. Slowly rotating stars 

Rotation affects the structure of tlie star and 
the spacetime, altering the mass-radius relation, 
the frequency of an orbit of given radius, and the 
radius Rms of the innermost stable orbit. 

The parameter that characterizes the impor- 
tance of rotational effects is the dimensionless 
quantity j = c J/ GM^ , where J and M are the 
angular momentum and gravitational mass of the 
star. The value of j that corresponds to a given 
observed spin frequency depends on the neutron 
star mass and equation of state, and is typically 
higher for lower masses and stiffer equations of 
state. For the spin frequencies ^ 300 Hz inferred 
in the kilohertz QPO sources, j ^ 0.1-0.3. For 
such small values of j, a first-order treatment is 
adequate. To this order, analytical expressions 
are available for the relevant quantities and one 
can prove the existence of upper bounds on the 
mass and radius [10]. However, the bounds must 
be computed numerically. 

To first order in j, the orbital frequency (mea- 
sured at infinity) of gas in a prograde Keplerian 
orbit at a given Boyer-Lindquist radius r is 

i^Kir, M,j) « [1 - jiGM/rc^f^^] vUr. M) (6) 

and the radius of the innermost stable orbit is 

i?ms(M,j)«[l-j(2/3)^/']i?L(M), (7) 

where z/^ and ^^"^ the Keplerian frequency 
and radius of the innermost stable orbit for a 
nonrotating star. Hence, to first order in j, the 
frequency of the prograde orbit at i?„is around a 
star of given mass M and dimensionless angular 
momentum j is (see [4] ) 

i^K,m. « 2210 (1 + 0.75j) (Mq/M) Hz . (8) 

As a result, for slowly rotating stars the upper 
bounds on the mass become 

M„,ax « [1 + 0.75jVspin)]MO,, (9) 

and 

i?max « [1 + 0.20j(j.spin)]<ax , (W) 

where j(vspin) is the value of j for the observed 
stellar spin rate at the maximum allowed mass for 
the equation of state being considered and M^^^^ 



and i?max 9.re the maximum allowed mass and 
radius for a nonrotating star (see above). Ex- 
pressions (^ and ( |lO|) show that the bounds are 
always greater for a slowly rotating star than for 
a nonrotating star, regardless of the equation of 
state assumed. 

Figure Ic illustrates the effects of stellar ro- 
tation on the region of the radius-mass plane 
allowed for spin rates ~ 300 Hz , like those 
inferred for the kilohertz QPO sources, and 
^QP02 — 1220 Hz , the frequency of the highest- 
frequency QPO so far observed in 4U 1636—536, 
which is also the highest-frequency QPO so far 
observed in any source. Our calculations show 
that the mass of the neutron star in 1636—536 
must be less than ^ 2.2 Mq and its radius must 
he less than ^ 17 km. As just explained, the 
precise upper bounds depend on the equation of 
state assumed. For further details, see [10]. 

2.3. Rapidly rotating stars 

If the stellar spin is ^ 500 Hz or higher, spin 
affects the structure of the star as well as the 
exterior spacetime. The exterior spacetime of 
such a rapidly rotating star differs substantially 
from the Kerr spacetime and must be computed 
numerically for each assumed equation of state. 
Derivation of bounds on the mass and radius of 
a given star for an assumed equation of state re- 
quires construction of a sequence of stellar models 
and spacetimes for different masses using the as- 
sumed equation of state, with Vg-pin as measured 
at infinity held fixed. The maximum and min- 
imum possible masses and radii allowed by the 
observed QPO frequency can then be determined. 

Such computations have been carried out by 
Miller, Lamb, & Cook [9]. They find that ff the 
neutron star is spinning rapidly, the constraints 
on the equation of state are tightened dramati- 
cally. For instance, if the spin frequency of the 
neutron star in 4U 1636—536 is ~ 580 Hz , the 
frequency of the single brightness oscillation ob- 
served during X-ray bursts from this source, then 
the tensor-interaction equation of state of Pand- 
haripande and Smith [13] is ruled out by the 
kilohertz QPO frequencies already observed from 
this source. They also find that observation of 
a 1500 Hz orbital frequency would constrain the 
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mass and radius of the neutron star to be less 
than ~ 1.7 Mq and ^ 13 km, ruling out several 
equations of state that are currently astrophysi- 
cally viable, regardless of the star's spin rate. 

2.4. Effects of the radial radiation force 

The luminosities of the Z sources are typically 
'--^ 0.5-1 Lb, where L^; is the Eddington luminos- 
ity. Hence, in the Z sources the outward acceler- 
ation caused by the radial component of the radi- 
ation force can be a substantial fraction of the in- 
ward acceleration caused by gravity. The radially 
outward component of the radiation force reduces 
the orbital frequency at a given radius. For exam- 
ple, if the star is spherical and nonrotating and 
emits radiation uniformly and isotropically from 
its entire surface, the orbital frequency (measured 
at infinity) of a test particle at Boyer-Lindquist 
radius r is given by 



ML) 

MO) 
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(1 -3GM/rc2)i/2 L 



{l-2GM/rc^) Le 



1/2 



(11) 



where L is the luminosity of the star measured at 
infinity and i'k(O) is the Keplerian frequency in 
the absence of radiation forces. Thus, the Boyer- 
Lindquist radius of a circular orbit with a given 
frequency is smaller in the presence of the radial 
radiation force and the constraints on the mass 
and radius of the star are therefore strengthened. 

For the atoll sources, which have luminosities 
< 0.1 Le, the change in the Keplerian frequency 
is at most 5%. For the Z sources, on the other 
hand, which have luminosities ~ Lb, the change 
may be much larger, although the change in the 
sonic-point Keplerian frequency may be smaller 
than would be suggested by a naive application 
of equation (JTl]), if a substantial fraction of the 
radiation produced near the star is scattered out 
of the disk plane before it reaches the sonic point. 

3. INNERMOST STABLE CIRCULAR 
ORBIT 

Establishing that an observed QPO frequency 
is the orbital frequency of the innermost stable 
circular orbit in an X-ray source would be an 
important step forward in our understanding of 
strong-field gravity and the properties of dense 



matter, because it would (1) confirm one of the 
key predictions of general relativity in the strong- 
field regime and (2) fix the mass of the neutron 
star in that source, for each assumed equation of 
state. 

Given the fundamental significance of the in- 
nermost stable orbit, it is very important to es- 
tablish what would constitute strong, rather than 
merely suggestive, evidence that the innermost 
stable orbit has been detected. Clear signatures 
of the innermost stable orbit include the following 
(see [10] for more details): 

(1) Probably the most convincing signature 
would be a fairly coherent, kilohertz QPO with 
a frequency that reproducibly increases steeply 
with increasing accretion rate but then becomes 
constant and remains nearly constant as the ac- 
cretion rate increases further. The constant fre- 
quency should always be the same in a given 
source. As shown in Figure Id, this behavior 
emerges naturally from general relativistic calcu- 
lations of the gas dynamics and radiation trans- 
port in the sonic-point model [5,10]. 

(2) A possible signature of the innermost sta- 
ble orbit would be a simultaneous sharp drop in 
the amplitudes of both QPOs in a kilohertz QPO 
pair, or a sharp drop in the amplitude of the 
lower-frequency QPO, at a frequency that is al- 
ways the same in a given source. 

(3) A less likely but possible signature of the 
innermost stable orbit would be a steep and re- 
producible drop in the coherence of both QPOs in 
a pair (or in the coherence of the higher-frequency 
QPO, if the lower- frequency QPO is not visible) 
at a certain critical frequency (the frequency of 
the innermost stable orbit), as the frequencies of 
both QPOs increase steadily with increasing ac- 
cretion rate. The critical frequency should always 
be the same in a given source. 

Several authors have recently suggested that in- 
nermost stable orbits have already been observed. 
Zhang et al. [22] suggested that the similarity 
of the highest QPO frequencies seen so far in- 
dicates that innermost stable circular orbits are 
being detected. Kaaret, Ford, & Chen [3]3 sug- 
gested that the roughly constant frequencies of 
the 800-900 Hz QPOs seen in 4U 1608-52 [2] and 
4U 1636—536 [21] were generated by the beat of 
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the spin frequency against the frequency of the 
innermost stable circular orbits in these sources. 
This would imply that the neutron stars in all 
the kilohertz QPO sources have masses close to 
2.0 Mq. However, no strong signatures of the in- 
nermost stable circular orbit have so far been seen 
in any of these sources. 

Indeed, more recent observations of both 
4U 1608-52 [8] and 4U 1636-536 [19] are incon- 
sistent with the suggestion that the QPO frequen- 
cies seen initially are related to the frequencies 
of the innermost orbits in these sources. The 
1171 Hz QPO frequency assumed by Kaaret et 
al. [3] to be the frequency of the innermost sta- 
ble orbit in 4U 1636—536 and used by them to 
determine the mass of the star was later seen to 
increase to 1193 Hz and then to 1220 Hz ([19]; 
W. Zhang, personal communication) . There is as 
yet no evidence for a maximum QPO frequency 
in 4U 1636—536 and hence there is no basis for 
the suggestion that an innermost stable orbit has 
been seen in this source. 

A recent analysis of 4U 1608—52 data by 
Mendez et al. [8] shows that this source has 
twin kilohertz QPOs that vary with countrate 
just like the other atoll sources. There is as 
yet no evidence for a maximum QPO frequency 
in 4U 1608—52 and hence there is no basis for 
the suggestion that an innermost stable orbit has 
been seen in this source, either. 

Even though no strong evidence for the effects 
of the innermost stable circular orbit has yet been 
seen in any kilohertz QPO data, there is reason 
for optimism, because the highest QPO frequen- 
cies observed so far are only 100-200 Hz below 
the expected frequencies of the innermost stable 
orbit around these neutron stars. Given the rapid 
pace of discoveries in this field, the prospects for 
obtaining clear evidence of an innermost stable 
circular orbit in the future appear excellent. 
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